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ABSTRACT 


шонстошсстс developed an analytical model of a repair- 
able item inventory system. The system consisted of a 
depot that repaired failed units according to some repair 
policy and stocked ready-for-issue units in support of a 
Zinite number Of customers. A least-cost repair policy 
and stock level was determined by use of a computer program 


taich is included as an appendix. 
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I. BACKGROUND 


This thesis studied the problem of inventory management 
of repair parts that have high unit cost and few customers. 

The importance of a model such as this lies in the fact 
that, although these parts represent a very low percentage 
of items managed by the Navy, they account for a large 
percentage of the dollar investment in inventory. 

Because of the high dollar value, it costs less to repair 
failed units than to buy new ones. Thus, it pays to establish 
а repair facility to renew failed items. 

This thesis attempted to model this repair process to 
determine how this type cf system might operate at minimum 


cost. 
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11. INTRODUCTION 


A repairable item inventory system consists of three 
ешешетсек users or customers, repair facilities, and 
stocking activities, 

When a part fails, the customer forwards that unit to 
the repair facility and requests a new unit from the stocking 
activity. Upon receipt of the failed unit, the repair 
facility refurbishes it. The number requiring repair at any 
one time is called the "repairable stock". The third element, 
Mets tbockine activity, stores the parts received from the 
repair facility and issues them on request from the end user. 
The number of ready-for-issue parts in the stocking activity 
is called the safety stock. 

This is the repairable item inventory system. One 
should note that the system as defined is closed. That is, 
every failed unit received by the repair facility can be 
repaired. Thus, after the inventory system is in operation, 
there is no reason to order new units to replace those 
already in the system. 

In order to operate the repairable item inventory system 
optimally, two basic questions must be answered. 

1. How many failed units should be in the repair 
Secility before one should repair and how many should one 


repair when the decision to repair is made. 





2. How much safety stock should be carried by this 
stocking activity so that orders from end-users might be 
promptly and economically filled. 

When one has made these two decisions one can compute 
a cost per cycle to operate under these conditions. The 
goal of course is to make these decisions so that the system 
operates at minimum cost. 

The answer to the first question, when and how much to 
repair, was determined using a technique developed by Howard 
in his Dynamic Programming and Markov Processes. It is, 

Ша наср, (пе purpose of this thesis, to look into the 
Bcusubsüsbty of utilizing such a method on this problem. 

The answer to the second question, how much to stock, 
can be found only after solving the repair problem. 

Since the problem is naturally divided into answering 
question (1) and then question (2), Section III will treat 
the problem in this sequence. In Section IV, there will be 
a sample solution and a discussion of the computer program. 
ama ly. cne methods used willbe criticized, and 


conclusions will be made in Section V. 
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A. GENERAL 
A model such as the one that was described in the 
introduction requires that a decision about the repair policy 
be made at specified time periods. The point in time at 
Which this decision is made is called the review time. The 
time between reviews is held constant and because of this 
the model is called a periodic review model. The length of 
time between reviews is defined as the average time needed 
to repair а failed unit. This is called the repair cycle. 
(ере по the length of the repair cycle in this way, the 
model with its associated assumpticns will better fit the 
real-world situation. This will be seen more clearly later. 
Graphically, the repairable item inventory system operates 


in the following manner during a repair cycle. 
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Customers Repair Depot 
Figure 1, 


REPATRABLE ITEM INVENTORY SYSTEM 
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Tn during the repair cycle, customers! repairable 
units fail and set in motion the dynamics of Figure 1. In 
this illustration there are N customers and M units of safety 
stock. 

Gemsider the effecte of this M and N upon the quantity 
of repairable stock. The maximum number that can be in 
meed of repair is N+M, while the minimum number is zero, 
гето за repair implies that each customer has a working 
unit and that there are M units іп a ready-for-issue con- 
dition in stock. At the other end of the range N+M units 
in repele implies that each customer lacks this parteandethat 
Mere ате zero units of safety stock available. This is the 
worst condition in which the system is allowed to function. 
Ше ount of repairable stocka a number between zero and 
N+M is defined as the "state" of the system. Thus, there 
are N+M states. 

ШИСИ ероод от courses would not need to repair all 
failed items on hand. The alternatives available are to 
repair from zero up to the total amount awaiting repair. 

For example, if there were 5 units of repairable stock, the 
e-peurcould Make Sax decisions: repair O, repair 1, ...; 
pe. Define this decision as the repair policy for 
that state. Define a set of decisions, one for every state, 
ав Бие "repair policy". 

The preceding has been a verbal description of a 
repairable item inventory system. The next few sections 


will translate this into mathematical terms. 


ШІ: = 





i. "Рев Ргобар 11 Бу Matrix 
Let S. be the amount of repairable stock at the 


th 
beginning of the n repair cycle. 


R be the amount of repairable stock actually 


repaired during the ntn repair cycle. 
Dr be the amount demanded during the neh repair 
Cycle. 
Thus S = – - Е == I: This describes a Markov 


enn. The distinctive element in a Markov chain’ is that 
the probability law of the future of a process, once it is 
in a given state, depends only on the state and not on how 
the process arrived in that state. 

Using the recursive definition of S above, a matrix 
can be formed that gives the probability of having any 
amount of repairable stock on hand, given the amount at the 
beginning of the previous repair cycle. This movement from 
one state to another is called a transition. 

Thus, if the present state is i, and the alternative 
chesen is to repair k, then the next state of the system 
will be j, if amdoniy if j+k-i items fail and hence are 
demanded during the cycle. 

Assumptions can be made concerning ЕЕ number and rate 
of these demands. A typical one might be that the demands 
are distributed Poisson with mean Du This assumption, 
although not necessary to the model, will be used in the 


remainder of the thesis. 
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If p, (1,3) denotes the probability that the current 
state is 1, k items are repaired and the next states is j 


then p (319) is given by: 


0 LE i-k< 0O 
JFK-1É0 
ү ОК | 
(1) py (1,3) 2 е A | dot aN-i-M & i AM 
| (j+k-i)! Stcherwise asi- O,..:., NIM 


К=0,....1 
J3=0,...,NM 


The restrictions upon p, (1,3) when p, (1,J)-0 in the 
formula above can be explained in the following way: 

а. i-k40. The probability p(.) equals zero when 
one wants to repair more than is in repairable stock. The 
decision is not a valid one. 

D. 110 This situation Can occur only if cher 
is negative demand, which is clearly impossible. 

с. jtk-i>N. This can occur only if there are more 
demands than customers can generate. 

de Шше-ізіеі М опа I SM. 

If iaM then the number of demands in the cycle can 
be no greater than N-(i-M) since i-M customers will not have 
items. Thus those cases in which j+k-i exceed N-i+M must be 
assigned zero probability. 

Suppose A is a vector that was defined earlier as 
шие repair policy, that is, a repair decision for en 
possible state. The transpose of the (N+M+1.) vector A 


Sen pe denoted by: 
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Сш. (a。， а))...» 24» .. PP 


where а, is theramount to be repaired, if the 
: $ „th 
Begeır cycle starts in the i. state. 
For any given repair policy, a probability matrix 


describing this passage from state to state can be defined 


= о М+М 
0 | esce A6 
то вар Ж шш. 
(2) ОШО pM. M СС 
pes . [7 
E тш л 


Wnere the P's denote those probabilities 
which are positive 


Because of the number of alternative repair policies 
available for a given starting state, there are (N+M+1)! 
possible probability state matrices. By defining the 
matrices in this way, three additional assumptions are 
made: 1) even after a decision to repair items in the 
repairable stock is made, those items will not be available 
until the next period. In other words, demands can only 
be filled from those items that are ready for use at the 
beginning of the repair period. 2) Parts that have failed 
during a repair period are delivered to the repair гасайа су 
prior to the end of the period regardless of when they fail. 
3) Shipments from the repair facility are received by the 


customers at the end of the repair cycle. == 
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2 нотрое се Етвойте ву 


To use the techniques discussed in this thesis, the 
states of she probability matrices described in the previous 
section must not just be Markov chains, but ergodic Markov 
chains. Taat is, the one-step probability transition 
Matrices which result from repair policies must be completely 
ergodic. By this is meant that after a large number of 
transitions, the state of the system is independent of the 
mal state. For a matrix P(A) оре харао а о санасу: 
it is sufficient to show that it has one recurrent chain. 

Consider first a repair policy in which the alter- 
native is "do not repair" regardless of the state of the 


system. The one-step probability matrix would be as follows: 


О т 2 N N+i N+M 
O PP PereP Oe... Ọ 
O o Oo О 
al О Po Der: "PS Бе о ыы то 
О 


210 «OL. Pape pp 


Figure 2, 


Repair Policy: "Do Not Repair" 


The symbol Ба denotes those probabilities which are 


positive under the repair decision "do not revair." From 





the figure one can see that eventually the system will end 
Up in state N+M, independent of the initial state, and 
remain there. 

To get & better feel for this property in conjunc- 
tion with & repairabie item inventory system, let the matrix 
shown in Figure 2. change slightly. Consider the probability 
transition matrix which results from changing just one of 


the alternatives: repair k items if the state is i. 


0 1 N+M 
О E Po ° e Po О Бо 
ni O ша : e Y е ES о O е О 
Іі: ES E Р, 2 | 
N+M E е е е е ° е е Po 
Figure 3. 
Repair Policy: "Repair k if in states i; 


otherwise do not repair any" 


The decision to repair k items shifts the associated 
МӨ ОГ the matrix k units to the left. Thus, a decision 
to repair is equivalent to allowing the system to shift 
back to some lower state with positive probability. This 


must be kept in mind when considering the next question. 
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ss there any repair policy that would make this 
probability transition matrix non-ergodic? Consider the 
following repair policy: If in states O through N+M-1, 


repair all; if in state N+M, repair none. 


О 1 2 N N+M 
О Po Ро PS . . ° Зе О. О 
dl Pz P4 P4 . ° . o P4 Q o Q 
2 Po Р 2 Р 2 » • > • Ро О ? О 
КОКЕ РЕ - = Ра RO 
m m m 
+ : : 
М+1 Patl nil mi Я ана 2 
N4M-1| - О 
N+M| O 0 О е : M е О P 
в nim. 
Figure 4. 


Non-ergodic Repair Policy 


im this example, if, after a long period of time 
one found himself in state N+M, one would know he had started 
in ЕМ. Thus, this repair policy produces a probability 
transition matrix that is non-ergodic. 

A non-ergodic repair policy could arise in three 
different instances in a problem of this nature. Each of 
these will be considered separately along with the method 


for handling that situation. 
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а. А non-ergodic repair policy might exist in the 
EC Of possible repair policies, but might not be utilized 
in any way while determining the optimal repair policy. 

The method developed in the following sections requires that 
only the repair policies considered be ergodic, not all 
possible repair @olicies. 

b. A non-ergodic repair policy could exist and could 
be in the set of repair policies considered in deriving an 
optimal policy that is completely ergodic. In this case, 
there is a method developed by Howard in his Dynamic Pro- 
gramming and Markov Processes that can solve this type of 
problem. 

con Einallys  non=sersodie repair policy might pE 
oea. But, ІГ а non=-ergodic policy were optimal, one 
would actually have two separate repairable item inventory 
problems. Thus, they could be separated and treated as such. 

From the example that had to be used to show that a 
non-ergodic repair policy might exist, one can see that it 
sen ınusVal case. So, with the solutions offered for 
these cases, no generality will be lost in considering all 
Meparr policies completely ergodic. 

3. Costs 

Шр until nom the thesis bhas described the repairable 
item inventory problem as a series of transitions, going from 
one amount of repairable stock to another, depending upon 


meme pair policy. But, from this, nothing can be said 
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about the relative merit of one repair policy over another. 
Inis is a function of costs. 

For instance, safety stock serves as an alternative 
to repairing a failed unit. Having an item available in 
safety stock allows one to delay repair of the failed unit 
and still offer immediate delivery. Since both of these 
alternatives offer delivery at once, the choice of one 
mMeuweed Of Operation over the other is dependent upon costs. 

Mére is also the alternative between immediate 
пей and delayed delivery. But, here again, one can 
attach a cost for immediate delivery and delayed delivery. 
hus. there is a set of trade-offs and each has an associated 
Come, It must only be determined which repair policy and 
stock level will allow this system to operate at minimum 
cost. 

For this repairable item model, charges will be made 
at the end of each repair cycle. At that time one will 
determine what has happened during the previous cycle and 
what is on hand at the moment and then make appropriate 
changes. The following costs will be considered while doing 
this: 

a. stock Level Costs 

b. Set-up and Repair Costs 

c. Additional Holding Costs 

d. Backorder Costs 

а. Stock Level Costs 

This cost ineludes. the rental or depreciation 


of a warehouse. The maintenance, taxes, and insurance are 


ШУ 





also considered in this charge. These costs are considered 
fixed after the size of the operation has been determined. 

It is assumed, in this model, that if it costs 
$K in fixed costs for a stocking policy of M units, it will 
cost $2K to have a stocking policy of 2M. 

But, after the stocking policy has been set, 
the cost, be it $K or $2K, will be charged no matter how much 
safety stock is actually on hand at the end of the repair 
eye le. 

The stock level M will initially be assumed to 
be weavenmend the optimal repair policy will bewdetermined 
for that given level. In this instance, the stock level 
cost is constant and independent of repair policy. So it 
Mecammot be used in this computation, 

Later im the thesis the optimal choice of M will 
овна red; and, at that time, the stock level cost willl 
enter the calculations. 

ре Additional Holding Cost 

At times the depot might have an amount in 
excess of the fixed stock level M on hand. This comes about 
when the amount demanded exceeds the safety stock. An 
additional amount is charged for each unit in the depot in 
excess of this safety stock to cover the cost of such things 
as additional warehouse space and overtime required to 
Bemdle these items. 

So, if S is the number in repair at the end 


o Ebe period; h is the holding cost per unit 


20 


Then 


(3) ia 


h(S-M) S7M 


с) Set- up amare pacos Es 

When the decision is made to repair some failed 
units, two costs are incurred: а set-up cost and a repair 
соз. 

The set-up cost is a fixed cost that is levied 
every time the decision is made to repair and is independent 
ОТ the quantity repaired. It is considered the cost to - 
tool up. 

The second eost is the charge that is made for 
pepairing each individual unit. This cost includes the labor 


Ama materials required to repair it. 


So let C = the set-up cost 

r= per оласи разг сос 

х = quantity repaired 
Тпеп 

Cr DX KEO 
(4) Cost = 
rep 
O otherwise 


а. Backorder Costs 
THis is a cost charged when the quantity demanded 
exceeds available safety stock. It is a penalty for loss 
of goodwill or sales. In the Navy's case, it might be a 
penalty levied for lack of readiness. 
In this model the penalty cost can only be levied 


if repairable stock exceeds M. If it is less than M there 


will be at least one unit of safety stock available. 
Eu m 


o 





So if the number in repair at the end of the 
repair period is denoted by S, and P is the backorder cost 


per unit, then 


О S£ M 
(5) Совт аск = 
P(S-M) S 7M 
The total variable cost per transition would 
then be the sum of the repair cost, the holding cost, and 


the backorder cost. 


В. DETERMINATION OF THE OPTIMAL REPAIR POLICY FOR A GIVEN 
STOCK LEVEL 
With the model, the transition probability matrix, and 
the costs defined as above, the minimum cost repair policy 
Noria. giverstock level can be developed. 
1. The Problem as a Markov Chain 
Let Б(п) be a vector wnose ith element is the 
probability that the state of the system will be in i after 
nerepair cycles. P(A), as has been defined before, is the 
one-step transition probability matrix for a given repair 
policy A 
If one knows the state vector b(0) and the Matrix 


Р(А) then the state of the system after one repair cycle 15 
(6) b(1)=5(0)P(A) 


If one wishes to know the state vector b(n+1) with 


b(0) known then 


22 
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(7) b(n+1) = b(n) P(A) 
5(n-1)[ 520] ^ = 8-2) [909] ? - 
P 


5(о) | Р(А) | 


(8) 


Now, if P(A) is ergodic, it can be shown that steady 
state probabilities exist and can be determined by solving 


the following equations: 
(9) b = bP(A) 


(10) 1'b = 1 
where the elements of b are the steady-state 
probabilities of having a given quantity in repair. 
suppose then that repair policies vector A consists 


or ай BODEN and suppose further that there are 


ЕОМ mern 


Pedros Of repairable stock, The expected cost for the 
period would be 
e ) К (3) 5 150 Ра, (5522,4. (isd) 
мһеге р (aj) is the probability of endins In 
ai state j, having started in 1, using 


the repair policy, repair а. 


С. (i,j) ie the cost aeeeciated with p (5) 
3 а. 
3 й 


Probabilities Pa (i,j) and costs C, (i,j) are known 
1 1 
апа, so, expected cost K (i) can be easily computed. 
al 
Since the expected cost of being in each state is 


known and the long-run probabilities of being in each state 


are known (Equations (9) and(10)), the long-run expected 





cost per period can be determined and would be: 


(12) K, = b Ko, (1) 
159 2 Б 

"ne could look at each possible repair policy, A, 
and compute this К Then the problem of determining a 
minimum cost policy for a fixed stock level M could be solved 
Dy cemparing each cost, = and choosing the minimum value. 
Bom ver, this would involve determining a cost for a large 
number of repair policies. 

For example, for three customers and a two unit 
Боса level, the number of repair policies that must be 
compared is 6! or 720. Thus, even for a very small size 
problem, the number of repair policies to be examined is 
апте large. For a larger problem the number of repair 
policies becomes so large as to prohibit even the use of a 
Gempuver to solve the problem in the manner described above. 

Thus the method is feasible but computationally not 
Pmaetical. The equations developed in this section, though, 
Will be used later. 

An alternative procedure for solving the problem 
imsasmenner which requires less computation is necessary. 

2. The Problem as a DRE Program 

To gain an understanding of this alternative approach 

пиле вез. to look initlally at a finivte time horizon. 


Consider a repairable item system over the next 


L-periods. Let there be n cycles remaining until the time 
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interval terminates. Define a, (n) as the repair policy 
that will be utilized if the system has i units in repair 
ema n cycles remaining. When i is specified for each n, 
the repair policy over the life of the system will be 
specified. That repair policy which minimizes the total 
expected cost over the entire time period would be optimal. 
Define V: (n) as the minimum total expected cost of 
maintaining the system with n cycles remaining, having started 
with i in repair (assuming an optimal policy has been followed 
us to this point). 
Vi (n) can be computed by first looking at the minimum 
costrpolicy if the system had one cycle remaining. If at 


the start of this period, one had i in repair, then 


НАЕМ 、 | 
(18) М (1) = Му р,(1,2)с,(1,2) 
к 3-0 


With V,(1) known one could look at DENS 


N+M (prob. of Cost of 
(1%) а > Transition || Transition 
i GEG To state j to Staten 
repairing К 

Minimum expected 

cost for week 1 

if in state j 

3 N+M QUA des 

(15) v(e) = um "E (1,2) [e3 0] 


k j=0 


In this manner one could work backward and: determine 


рн орсішал policy for each cycle. 


(Ze 





Perecuncwve relation arises out of this approach. 


N+M | 
(16) ve) em Y se [e 0:2) + v 062)] 
j=0 


If termination of this repairable item system 
were imminent, this approach would be valid and efficient. 
If one did not know the number of cycles until termination, 
but only that it was a long time into the future, this method 
seems less effective. 

Consider, then, what happens to V: (n), defined 
above, when n gets large. To accomplish this goal, 
generating functions will be used. 

3. Asymtotic Form of Recursive Relation 

The generating function of an arbitrary function 

defined on the non-negative integers is defined by 


ог 
n 
(17) Sz) mn) а 
п+0 
Some general properties of generating functions 


are the following: 


General Properties of Generating Functions 


Generating Function 


iz 
ЛЕНЕ 


z JE (z)-£(0)] 





ТАВИН Е 
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During the following development, generating functions 
of vectors will be used. A vector generating function is 
Simply a vector whose elements are generating functions. 

To develop the asymptotic form of V, (2) (Equation 16) 
through generating functions, one must also develop the 
asymptotic expression for b(n) (Equation 7). 

Transforming Equation Y using the above table 


Produces: 
а = АМЕ = 
(18) 27 [5(2) - 5(0)] = 5(2) Р (А) 
With some manipulation, the above equation becomes: 


7 = = 二 SR E 
(19) b(z) = b(0) (I - zP(A)) 
OttintGerest in this equation, for use a little later 
ІШ the section, is the inverse transform of the following 


matrix from Equation (19) 
= = == Al 
(20) Н(а) = (1 - = Р(А)) 


re 
Consider (Ii-zP(A)) . By the adjoint method of 


emaluatine an inverse, one can write: 


(21) H(z) = Adjoing (1-2) 
(I - 2Р(А)) 


sance z=l is a characteristic root of every Markov 


Metrix, {Т-®Р(А)) can be written as: 


(22) (1-2) оа чш 7 Ыга ) 





Thus by partial fraction expansion of Equation 21 


one can show that H(z) = S + G(z) 
pz 


where G(z) is & rational function of z 
Since from Table I, the inverse transform for 


S/(1-z) is 5, the inverse transform of H(z) would be 
(24) H(n) = S + G(n) 


H(n) is simple (ғ(А) |" as can be seen by 
Equation (8) and Equation (19). Since P(A) is a completely 
ergodic matrix, it must be true that G(n) > O авп->0 
and S is the steady-state probability matrix. 

зо. each rom of S is the row vector b. 

with this information at hand, an equation for 
V(n) een be developed. Suppose now that the repair policy 
has been given. The subscript, A, describing this polrcy 
will be dropped from the development since it adds no 
information. Multiplying through by p(i,j) Equation (16) 


becomes: 
М+М N+M 
(25) кар = a pie (2590 £ = BIT, UE (n) 
г ae 359 Б 


1=1,...,NM 
Combining Equation (25) with the notation of Equation (11) 
this becomes 
N+M | 
(26) VA) И A 
E ј=о Ј 


Switching to matrix form this can be written: 
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em) V(n+1) = K + PV(n) 


where КЬ is a column vector of the K(i)'s. 
Through the use of transforms in Table I and the 
fact that the transform of f(n) = 1 is 1/(1-z), the matrix 


of generating functions becomes 


(28) 277 [¥(2) - (0) | = 1/(1-2) + P V(z) 


With a little manipulation this becomes 
(29) V(z) «4 2-2 P) Ik, € (1-22) 1v(0) 


Substituting Equation (24) into Equation (29) and collecting 
terms, one gets: 
(30) T(z) => „SK +Z G(z)e i SV(0)4G(z)V(0) 
(l-z) b 1c l-z 

Since the inverse transform of z/(1-2)° is f(n)=n, 
mer the first term ison SR. By using partial fraction 
expansion and dropping terms that tend to zero as n becomes 
large (the asymptotic form is only of interest), the inverse 
transform of the second term is G(1)Kp. The third term 
bas an inverse transform of SV(0), since the inverse 
transform of 1/(1-z) is f(n) = 1. Finally, since the 
geometric terms of G(z) V(O) approach zero as n gets 
large, this term can be ignored. Thus the asymptotic 


inverse transform of V(z) is: 
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(31) V(n) = nSK, + G(1) К, + ОО) 


Since the matrix S is made up of N+M b vectors, 


one can write: 


(32) K = SK, 
where each element of K is the constant КА рег 
Equation (12). 


So V(n) becomes 
(33) V(n) = nk + G(1) К + 8 V(0) 


What has been derived here is an equation for long- 
run cost as a function of N. Looking at the equation in 


this manner, G(1)X, +sv(0) will be defined as V, the 


ексер. Ine asymptotic Iong-run cost equation finally 
becomes: 
(34) . V(n) 2 nK V 


4. Review and Collection of Information 

To quickly review, one initially wanted to find a 
method of determining the minimum cost repair policy. To 
do this, it was decided first to determine the minimum cost 
Bepeir policy given a certain stock level. 

Two approaches were used and combined: the Markov 
chain and Dynamic Programming approach. In the former, an 
dation for the cost per cycle can be obtained but it is 
shown that it is computationally impractical to solve a 


problem in this manner. In the latter approach, a simple 
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asymptotic form of total expected cost after n repair 
cycles is found. This leads to a simplified method of 
solution. 

Suppose the repair system is operating under a 
СЕО ролі су. Тет. by breaking the asymptotic form of 


the matrix V(n) into its elements one gets; 
(35) Vi(n) 9» nK * V, | i-0,...,N4M 
From equation (26) one has: 

T 


(36) Cc ET si p(i,3)V,(n-1) 1-0,...,N4M 
ј=0 


Now combining these two equations, 


N+M 
nK+V, = K(i) + ) p(i,J)V (n-1) 
(37) т 9-9 J 
N+M 
пК+/, = К(1)+ L с 
Ез 
N+M 
Шемечасе ) р(1,2) -1 
TO 
N+M 
(38) ЕЕ = ae) Lpi) JO NM 
=o 


Equation (38) produces N+M+1 simultaneous equations with N+M+2 

unknowns. That is, there are N+M+1 V;'s and the K as unknowns. 
But, it one of the E s are веске пал со zero, one 

can solve these equations. This causes all of the NE to 

be transformed by a constant amount, but does not affect the 

Wo snow this, let: 

(39) ШІ. = У. + а 


aL d 


в 


- — 
u 
= Е а. 
E 





Then 
N+M 
ки Rd) + p pen va 
(40) к. 
КК Фа = K(i) + È, pies) [v +a] 


N+M 
KV = K({i) + ) p(isd) Y, 
a 10 J 


Thus, scaling down the V's by setting one of them 
Gamal to zero does not affect the computation of the cost K 
of the system. 

The development of these formulas assumes an 
arbitrary repair policy and thus computes the cost per 


Hepear cycle for this policy. Equation (12) was an equally 


easy way to solve this identical problem. The advantage to 


the present way, though, is that it presents & method of 
Seertcing with an arbitrary repair policy and rapidly con- 
ПИШЕ "со the optimal repair policy that the previous 
method did not provide. 

5. Alsorithm 

Equation 16 showed that if we were in repair cycle 

Buche sleeast cost policy forzeyele n+tl in the ith state would 
pe 


N+M 
(41) MIN Е + е 2854) vn) | 


ај 


Since in the long run 


(42) V4(n) =пК+ у, 
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and 
N+M 

(43) 2 р (i,j) - 1 
j=0 


we can minimize cost over possible decisions in the Де state 
by choosing the decision: 
N+M = 
(шы) em [,, 0 +X, Pay (tod) Vy | 
Using Equation (44) in conjunction with Equation 
(38) produces an interactive scheme that converges to the 


optimal policy. This will be proved in the next section. 


Assuming this, though, the algorithm would operate as follows: 


ALGORITHM | 


Ше Far an arbitraryeset of decisions, or for decisions 
generated from Part B, solve the set of simultaneous 
equations (setting one of the V.'s = 0). 


e) 
(38) N+M 
K + Vi = K(i) + Y Diss) E: i=0,...,N+M 


ј=о 
В. Using the V_'s obtained in the above equations, 
compute over all T- and all а; '5 
(44) | М+М 52 
Кы, (1) + 2, ра; (1:4) ү, 
For each state choose that repair alternative Гог 


which the equation above is minimum. Go to part А. 
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Hach iteration will decrease the value of K. Upon 
reaching the optimal policy the same decision will recur. 
Thus, the system will have converged to the optimal repair 
policy and the variable cost for a given stock level. 

Now, it only remains to prove the point, and a 
|ui method of Solving a repair problem” of this nature 
Will have been found. 

6. Proof of Convergence 

Suppose we have solved the system of simultaneous 
equations using a set of repair decisions A. Suppose further 
that in using part B of the algorithm, a repair policy B 
is generated. If it can be shown that K Ky» convergence 
by this iterative scheme will be proven. 

Since policy B produces & lower cost than the policy 


Bonon part B of the algorithm, we know that 


N+M Қа А N+M WT 
(45) K(i) TUNES Pp(i: J) у. < K, (4) у p„(is3)V; 
Uo j=0 
502222211 
et 
М+М А ММ Т 
(16) C =K (1) + L pA(1,3)V; 7 Y pg(,3)V, -Kg(i) 
Jes J=0 


Thus с. IS greater than or equal Тош о; 


From part A of the algorithm we know 


B | N4M d B 
(47) Ка + У, = Кра) + y pg (1,3) V4 
j=0 


= я 





(48) А A l N+M A 
p E KP) + X Palid) | 
j=0 


Subtracting the above two equations results in 


кеа МБ = KG "E rio n. 
Ld X uU) EE Le [ 2 (1,3), -p, (1,5) 


, 


= 
ee) 
М+М д М А 
= = » РАЈ), 十 y рҺ(8520У; 
1-0 3-0 
М+М В М+М в 
FE Pa A Први 30 
Шы. ј=о 
М+М А 
-c ве p (3,3) v Pv | 
WE 
j= 
Let 
B 
(50) A(K) =K" -K 
Bays v Ay 
1 1 
So the above equation by substitution becomes: 
N+M 
ү K) + V =C + | V 
(51) AG) * A(,) 7 6, * E 54,3) Á(V,) 
1-0 


Equation (51) is identical to Equation (38) except 


that the terms are differences rather than absolute quantities. 
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Thus by analogy the solution to Ax is 


(52) Мк) = y wc 


j-0 


18 

where bj is the 1% element of b from Equation 9. Since all 
the elements of b are greater than or equal to zero, and 

C > O by definition, ÄAk>o. SoxP is less than K^ and 
the algorithm converges. 

So for a givenstock level a method for determining 
the optimal repair policy and the variable cost in a few 
iterations has been found; its value increases as N+M+1, 
increases. One must now develop a method of determining 
ӘНЕ” optimal stock level. 

7. Computation of the Total Cost 

In paragraph III.A.3.a., it was pointed out that 
there is & cost associated with the stock level that was 
fixed while determining repair policy. This fixed cost 
(FC) is added to the variable cost (VC) for & total cost 
for a given stock level. 

TC = FC + VC 
There remains to determine the stock level with its 


repair policy that gives minimum cost. 


С. DETERMINATION OF AN OPTIMAL STOCK LEVEL 
In determining an optimal stock level, the range of 
stock level can be any number between zero and infinity. 


One cannot, therefore, simply compute the optimal repair 
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policy for each stock level and then choose that stock level 
that generates the minimum total cost per cycle. 

The total cost is a function of stock level and cannot 
be written without considerable computation. So, to examine 
jie Navure of the cost Tunction, it will help to analyze 
the fixed costs and variable costs associated with each M 
separetelyy Fixed cost is a linear function of M so, it 
Сащсез по particular problem. Variable cost though requires 
a little more analysis. 

1. Analysis of Variable Cost 

a Set-Up Cost. 

As the stock level increases there will be 
fewer occasions for which one will repair failed items while 
Meaamizing cost. Thus the set-up сов: реу сус бе асс 
Zero as stock level approaches infinity. 

Б. Backorder cost. 

The greater the stock level the less frequently 
the system will enter a backorder state. So, again, the 
Backorder cost per cycle will tend to approach zero as the 
stock level approaches infinity. | 

с. Repair Costs. 

The costmbe repair the failed items is indepen- 
dent of the stock level. In the long run, each item which 
fails will be repaired. Since the expected number of 
ет тег рехг cycle isA, the expected repair cost per cycle 


must be Atimes the cost of repairing the item. 


зи = 





Adding the three costs mentioned above, one can 
па щото, аз stock level increases to infinity the variable 
Secueper cycle approaches the average repair cost per cycle. 

2. Analysis of the Total Cost 


The graph of the total cost functions, can be plotted 


as is shown below: 


p Total Cost 


OI ARO Cost 








Variable Cost 


¡Mean Repair Cost 


PURO. 
GRAPH OF THE TOTAL COST FUNCTION 


As an example of this total cost function the 
following parameters were used to solve a problem of this 
nature: 


Number of customers: T 
Mean Cyclic Demand 

Set-Up Cost 

Repair Cost 

Backorder Cost 

Additional Carrying Cost 

Per Unit Fixed Carrying Cost 


| ДОО ГО О 
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Total Variable Total 





ock Leve 
Sit L 1 Ba Fixed Cost Total Cost 
1 14.24 1:00 15.24 
2 11.29 2 .00 13.27 
3 9.45 3.00 12.45 
д 8.52 1200 12.52 
5 7.90 5.00 12.90 
6 152 6.00 13.52 
7 oer 7.00 тот 
8 7.09 8.00 15207 
9 6.93 9.00 19.95 
10 6.82 10.00 16.82 
TABLE ТТ 


` Неге, the variable cost is approaching $6.00 while 
fixed cost is a linear function of stock level with 1 as the 
5 оре. 
The heuristic arguments presented here indicate that 
wen cotal costs function is convex. The importance of this 
mies in the fact that if there is a stock level let us say a, 


such that TC, is the total cost of this stock level and if 


тс TC 
a а+1 


апа 
TO, X. TC, 4 


then Ед is the minimum total cost repair system and there 15 


met Tc that has a lower cost in the range. 


a) 
o Sma <00 





ТУ. AN EXAMPLE 


In the example initially worked by the author, the number 
of customers was small and as a result, all possible repair 
policies as discussed in Section III.A were studied. But, 
even with all these alternatives, each state eventually 
converged to onesot two repair policies; repair none or 
repair all. This seemed reasonable since, with a set-up 
cost, if one repaired any, it would be the least-cost 
очну to repair all failed unite. 

So, in expanding the size of the program, all other 
possible repair policies were deleted. Poisson demands 
were also initially assumed in the thesis but as the mean 
number of demands per cycle increased the calculation of 
the Poisson probability became computationally unfeasible. 
To cireumvent the difficulty the following theorem was used: 

Theorem: A random variable K that is distributed 
Poisson with a mean, m, approaches the Normal distribution 
with mean, n, and variance ms as m gets large. 


1 -3y 
lim = fe d 
m- oo к! ү27{ ss 


k-m-3 


|| 
o 


For m reasonably large, this theorem was used to 
approximate the Poisson probabilities in the computations 
discussed in this section. 

Golden section search was utilized to fird the optimal 


stock level. Briefly this method of search operates in tae 
小 0 = 





foliowing шаппег: 


AA ри ии 


а, b e 


Let FI be the range of search. Let b and € be points 
picked in some efficient manner (golden section search). 
nen: 

I. compute costs for b and c. 

2. if the cost at bé the cost at с, search new interval 
a,d using the same method. 

EE пе the cost at b7 the cost st €, search new interval 
ЕР Busing the same method. 

A method such as this will allow one to get as close 
as one wants to the optimal or converge to the optimal, 
depending on the method of choosing the points b and C at 
each stage. 


A computer program has been written and can be found in 


the back of the thesis. The program utilized the following 


data: 
Number of customers de 
Mean cyclic demand ШЕ 
Set-up cost 20% 
Repair Cost 3. 
Backorder cost ES 
Additional carrying costs e 
Ber unit fixed carrying cost ar 
Stock level range 0-24 


Using this data the computer routine picked the stock 
levels in the sequence shown from left to right and 


generated the following rerair policies: 
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The program converged on the optimal policy of maintaining 
a stock level of 22. In addition it made the following 
пера г decision: if, at the beginning of the repair cycle, 
the number in repair is greater than or equal to 8, one 
should repair all the failed units. Otherwise, do not 
repair. | 

This program used 236K bits of computer space and 


6 min.30 second of computer time. 
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V. CRITIQUE AND CONCLUSIONS 


A. SUMMARY 

Howard's computation technique was used in determining 
the optimal repair policy and the variable cost for a given 
stock levei. The algorithm associated with this can be 
мола in Section III.B.5. The total cost for this stock 
iel wes then computed by adding its fixed cost to the 
variable cost. 

It was then argued that the cost of operating the repair 
system was a convex function of the stock level. Thus, 26s 
of the many unimodal search methods could be used. With 
this information at hand, the following algorithm could be 


used. 














Minimum Cost Algorithm 


IS les шо stock levels overta general ransezueme 
a unimodal search method. 


2. Compute the total cost(s) using the repair policy 
algorithm and compare. Select range associated with that 
stock level that has minimum total cost. 


3. Choose new stock level to compare with present 
optimum stock level. 


4, If there is no stock level better, stop. Other- 
wise, go to step 2. 


44 E 


ша 





BR CRITIQUE 

Although this algorithm did locate the optimal repair 
policy and stock level, it presented some computational 
problems. In addition to this, there is some question as 
to the usefulness of the results because of the assumptions 
that were made. 

1. Assumptions 

The model envisaged a cyclic repair system in which 
all transactions took place within that cycle. At the same 
time, probability distributions were assumed that were not 
consistent with this assumption. For instance, a Poisson 
distribution does not prevent a part from failing a couple 
ӘН Geeonds prior to the end of the cycle. But, it is 
assumed that, if tris happened, the failed item would be 
delivered to the repair"“activity prior to termination of 
the eycle. 

It was further assumed that all work started in a 
epale cycle would be completed in that cycle. This implied 
that everything required to repair the failed item would be 
on hand at all times. This particular assumption might be 
good or bad depending upon what is required to repair the 
item, 

It was also assumed that every failed item received 
GewidNbe repaired. This is very unrealistic. But, it is 
felt that this assumption could be relaxed by allowing for 
some attrition, although no suggestions are to be offered 


at this. time. 
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There are two properties of this model, though, that 
make it particularly useful. First of all, the model assumes 
that the repair periods are cyclic. This Ts the way most 
workloads are formed and work is accomplished today. Second!y, 
the model is independent of the demand distribution used. 
So, even some empirical demand distribution could be used 
as easily as the Poisson. 

2. Computational Problems 

The two major weaknesses of the model are program 
ze Ja. Computer run-time. In the example in Section IV, 
it was noted that 236K bits of computer space were required 
for a problem that contained only 75 customers and a range 
of possible stock levels of 0-24. 236K is considered a 
very large program, while (5 Customers ls a very small 
problem. The large storage space requirement is due to 
the matrix size required to hold all the state probabilities 
HO. the associated decisions, In the example in Section ТУ, 
the matrix was 100x100x2 which uses 80,000 bits of computer 
space by itself. A single computer could not handle a 
problem with 1000 customers which is not an unrealistic 
tp rable item problem. 

The second major weakness was the running time of 
the program. The small example took 63 minutes. This 
lengthy running time was due to the fact that a probability 
of demand must be computed for each state and each repair 
alternative. And then, having done this, one was required 
to solve a series of systems of simultaneous =quations to 


determine the optimal repair policy. 
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3. "Bansmeter Sensitivity 


The parameters or costs were not tested in the 
thesis. This would normally be done to see how repair 
policies and stock levels were affected by changes in the 
set-up cost, and the unit costs of carrying and repairing 
failed items. It was felt, thoushs that the computational 
problem mentioned in the previous section required a solution 


prier to testing parameters. 


С. ALTERNATIVES 

There are ways to avoid some of the problems. The 
suggested alternative ways given here deal with either 
reducing the number of states or reducing the number of 
alternatives. 

i. Reduction of Storage Area 

Information could be generated as required for 

sach iteration rather than storing it. “Био, by doing coise 
computation time would be increased significantly. For 
example, even in determining a repair policy for a given 
stock level, the entire state space and cost space would 
have to be computed at each iteration. And, as was pointed 
out earlier, computer running time was already quite high. 
Thus, there does not seem to be a practical means of reducing 
storage space utilized while using the same approach 


developed in the thesis. 
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2. Compression of State Space 
A study of the results of Section IV led to two 


possible alternatives for reducing state space. Both 
alternatives required the following facts to be noted from 
WE Ith in Seetion IV. 

Let x be the amount of repairable stock on hand; 
eG ebesthe smallest x such that the repair decision is: 
repair all failed items. Then, if, at the beginning of a 
Repair cycle, x is greater than b, the decision will be 
Co repair all failed units. 

Secondly, note that the "active" state space is 
quite small. By "active" is meant those states or repair 
decisions that actually change during cost reduction routines. 
Tne size of the active state space seems to be a function 
more of cyclic demand and stock size than of the number of 
customers. The following two alternatives are offered: 


a. Reduce the Size of the Problem: 


Let M = the number in stock 
Amean = the cyclic demand 
N = number of customers 


Define the state space, the number in repair, 
in the following way. Let the state space = 0O,1,...,AMEANHM,K 
where K is any possible state greater than M+Amean. The 
bulk of the states of a larger problem are combined into 
еШ Ее К, A major difficulty in solving this problem is the 
determination of a cost associated with the state K. 
Oe Rouch no specific solution will be offered in this papers 
it is felt that this can be solved Meine expected cost 


racher tuan actual cost. 
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b. Reduce the Number of Alternatives for Certain States 
Define M, Amean and N in the same way, use the 
following general rule. For the states O,...M+Amean, one 
can make two decisions: repair none or repair all. For 
all states K greater than M+Amean, repair all. This would 
reduce storage space while retaining the structure of the 


problemn as it is. 


D. EXTENSION-MULTI-ITEM SYSTEM 
A discussion of a repairable item inventory system 
would not be complete without at least a cursory look at a 
repairable item system that is multi-item. Since this was 
not the main thrust of this thesis, a solution to this 
problem will only be suggested. 
1. “No Resource Constramt 
With the absence of constraints on resources, the 
muiti-item problem reduces trivally to a set of individual 
Meowarrable-item problems. The minimum cost repair policy 
and stock level can be determined for each item individually. 
2. Resources Constraint 
The more common repairable item systems would have 
a combination of two constraints: a budget constraint and 
a workload cons traint. Suppose Minimizing whe probability 
of a stock-out is the measure of effectiveness. To work 
with this measure, stock-out cost must be deleted from the 
cost function and the problem must be reformulated in the 


following manner: 
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i 
MIN i a, P(S/0), 


L 
< 
ше cj (BED 


x E. ү l F 
= 
where L = number of customers 


ЦЕ! 
P(S/0), = probability for the i item being 
T out ог stock 


а, = werenting factor for item т 


C.(P,) = cost of the repair policy associated 
a with item i 


E, = work units required to repair itemi 


x, = total number of the S item to be 
1 repaired 


D = total budget constraint 


F = workload constraint 


with the time required to compute the optimal repair 
су and stock level even for one item, solving this) type 
ош Minimization problem is not amenable Фо solution at this 


time. 


Е. CONCLUSIONS 

@aree major conclusions сап be drawn from this paper. 

ШИШ е total cost of maintaining a part in a repairable 
Шеш system is a convex function of stock level. 

2. The optimal repair policy will be of the following 


Рогш: 





bet x be the number of failed items in repair. Пе 
b be the state of the smallest quantity in repair for which 
КШ Wail be repaired.) Imen the following decision rule will 
be used: 
If xb, do not repair 
x2b, repair all. 

Mis is similar to the (s,S) policy in inventory theory. 

3. The dynamic programming-Markov process approach 
to the repairable inventory system leads to a minimum cost 


solution. 
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